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We propose a model for the dynamics of a limit order book in a liquid market where buy and sell orders 
are submitted at high frequency. We derive a functional central limit theorem for the joint dynamics of the 
bid and ask queues and show that, when the frequency of order arrivals is large, the intraday dynamics of 
the limit order book may be approximated by a Markovian jump-diffusion process in the positive orthant, 
whose characteristics are explicitly described in terms of the statistical properties of the underlying order 
flow. This result allows to obtain tractable analytical approximations for various quantities of interest, 
such as the probability of a price increase or the distribution of the duration until the next price move, 
conditional on the state of the order book. Our results allow for a wide range of distributional assumptions 
and temporal dependence in the order flow and apply to a wide class of stochastic models proposed for order 
book dynamics, including models based on Poisson point processes, self-exciting point processes and models 
of the ACD-GARCH family. 
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1. Introduction 

An increasing proportion of financial transactions -in stocks, futures and other contracts- take 
place in electronic markets where participants may submit limit orders (for buying or selling), 
market orders and order cancelations which are then centralized in a limit order book and executed 
according to precise time and price priority rules. The limit order book represents, at each point in 
time, the outstanding orders which are awaiting execution: it consists in queues at different price 
levels where these orders are arranged according to time of arrival. A limit new buy (resp. sell) 
order of size x increases the size of the bid (resp. ask) queue by x. Market orders are executed 
against limit orders at the best available price: a market order decreases of size x the corresponding 
queue size by x. Limit orders placed at the best available price are executed against market orders. 

The availability of high-frequency data on limit order books has generated a lot of interest in 
statistical modeling of order book dynamics, motivated either by high-frequency trading applica- 
tions or simply a better understanding of intraday price dynamics (see Cont (2011) for a recent 
survey). The challenge here is to develop statistical models which capture salient features of the 
data while allowing for some analytical and computational tractability. 

Given the discrete nature of order submissions and precise priority rules for their execution, 
is quite natural to model a limit order book as a queueing system; early work in this direction 
dates back to Mendelson (1982). More recently, Cont, Stoikov and Talreja Cont et al. (2010b) have 
studied a Markovian queueing model of a limit order book, in which arrivals of market orders and 
limit orders at each price level are modeled as independent Poisson processes. Cont and do Larrard 
(2010) used this Markovian queueing approach to compute useful quantities (the distribution of 
the duration between price changes, the distribution and autocorrelation of price changes, and the 
probability of an upward move in the price, conditional on the state of the order book) and relate 
the volatility of the price with statistical properties of the order flow. 

However, the results obtained in such Markovian models rely on the fact that time intervals 
between orders are independent and exponentially distributed, orders are of the same size and that 
the order flow at the bid is independent from the order flow at the ask. Empirical studies on high- 
frequency data show these assumptions to be incorrect (Hasbrouck (2007), Bouchaud et al. (2002, 
2008), Andersen et al. (2010)). Figure 1 compares the quantiles of the duration between order book 
events for CitiGroup stock on June 26, 2008 to those of an exponential distribution with the same 
mean, showing that the empirical distribution of durations is far from being exponential. Figure 
9 shows the autocorrelation function of the inverse durations: the persistent positive value of this 
autocorrelation shows that durations may not be assumed to be independent. Finally, as shown 
in Figure 2 which displays the (positive or negative) changes in queue size induced by successive 
orders for CitiGroup shares, there is considerable heterogeneity in sizes and clustering in the timing 
of orders. 

Other, more complex, statistical models for order book dynamics have been developed to take 
these properties into account (see Section 2.3). However, only models based on Poisson point 
processes such as Cont et al. (2010b), Cont and de Larrard (2010) have offered so far the analyt- 
ical tractability necessary when it comes to studying quantities of interest such as durations or 
transition probabilities of the price, conditional on the state of the order book. It is therefore of 
interest to know whether the conclusions based on Markovian models are robust to a departure 
from these simplifying assumptions and, if not, how they must be modified in the presence of other 
distributional features and dependence in durations and order sizes. 

The goal of this work is to show that it is indeed possible to restore analytical tractability without 
imposing restrictive assumptions on the order arrival process, by exploiting the separation of time 
scales involved in the problem. The existence of widely different time scales, from milliseconds to 
minutes, makes it possible to obtain meaningful results from an asymptotic analysis of order book 
dynamics using a diffusion approximation of the limit order book. We argue that this diffusion 
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approximation provides relevant and computationally tractable approximations of the quantities 
of interest in liquid markets where order arrivals are frequent. 



J QQ Plot of Sample Data versus Distribution 
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Figure 1 Quantiles of inter-event durations compared with quantiles of an exponential distribution with the same 
mean (Citigroup, June 2008). The dotted hue represents the benchmark case where the observations 
are exponentially distributed, which is clearly not the case here. 




0.5 



1 1.5 
events 



2.5 
ii vf 



Figure 2 Number of shares per event for events affecting the ask. The stock is Citigroup on the 26th of June 2008 



As shown in Table 1, most applications involve the behavior of prices over time scales an order 
of magnitude larger than the typical inter-event duration: for example, in optimal trade execution 
the benchmark is the Volume weighted average price (VWAP) computed over a period which may 
range from 10 minutes to a day: over such time scales much of the microstructural details of the 
market are averaged out. Second, as noted in Table 2, in liquid equity markets the number of 
events affecting the state of the order book over such time scales is quite large, of the order of 
hundreds or thousands. The typical duration tl (resp. Tm) between limit orders (resp. market 
orders and cancelations) is typically 0.001 — 0.01 ^ 1 (in seconds). These observations show that it 
is relevant to consider heavy-traffic limits in which the rate of arrival of orders is large for studying 
the dynamics of order books in liquid markets. 
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Regime 


Time scale 


Issues 


Ultra-high 


~ 10^3-0.1 s 


Microstructure, 


frequency (UHF) 




Latency 


High 


~ 1 - 100 s 


Trade 


Frequency (HF) 




execution 


"Daily" 


~ 10^ - 10^ s 


Trading strategies, 
Option hedging 



Table 1 A hierarchy of time scales. 





Average no. of 


Price changes 




orders in 10s 


in 1 day 


Citigroup 


4469 


12499 


General Electric 


2356 


7862 


General Motors 


1275 


9016 



Table 2 Average number of orders in 10 seconds and number of price changes (June 26th, 2008). 

In this limit, the complex dynamics of the discrete queueing system is approximated by a simpler 
system with a continuous state space, which can be either described by a system of ordinary 
differential equations (in the 'fluid limit', where random fluctuations in queue size vanish) or 
a system of stochastic differential equations (in the 'diffusion limit' where random fluctuations 
dominate) (Iglchart and Whitt (1970), Harrison and Nguyen (1993), Wliitl (2UU2)). Intuitively, the 
fluid limit corresponds to the regime of law of large numbers, where random fluctuations average 
out and the limit is described by average queue size, whereas the diffusion limit corresponds to the 
regime of the central limit theorem, where fluctuations in queue size are asymptotically Gaussian. 
When order sizes or durations fail to have finite moments of first or second order, other scaling 
limits may intervene, involving Levy processes (see Whitt (2002)) or fractional Brownian motion 
Araman and Glynn (2011). As shown by Dai and Nguyen (1994), there are also cases where such a 
'heavy traffic limit' may fail to exist. The relevance of each of these asymptotic regimes is, of course, 
not a matter of 'taste' but an empirical question which depends on the behavior of high-frequency 
order flow in these markets. 

Using empirical data on US stocks, we argue that for most liquid stocks, while the rate of arrival 
of market orders and limit orders is large, the imbalance between limit orders, which increase queue 
size, and market orders and cancels, which decrease queue size, is an order of magnitude smaller: 
over, say, a 10 minute interval, one observes an imbalance ranging from 1 to 10 % of order flow. 
In other words, over a time scale of several minutes, a large number A'^ of events occur, but the 
bid/ask imbalance accumulating over the same interval is of order ^/N <^ N. In this regime, random 
fluctuations in queue sizes cannot be ignored and it is relevant to consider the diffusion limit of 
the limit order book. 

In this paper we study the behavior of a limit order book in this diffusion limit: we prove a 
functional central limit theorem for the joint dynamics of the bid and ask queues when the intensity 
of orders becomes large, and use it to derive an analytically tractable jump-diffusion approximation. 
More precisely, we show that under a wide range of assumptions, which are shown to be plausible 
for empirical data on liquid US stocks, the intraday dynamics of the limit order book behaves like 
as a planar Brownian motion in the interior of the positive orthant, and jumps to the interior of 
the orthant at each hitting time of the boundary. 
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This jump-diffusion approximation allows various quantities of interest to be computed analyti- 
cally: we obtain analytical expressions for various quantities such as the probability that the price 
will increase at the next price change, and the distribution of the duration between price changes, 
conditional on the state of the order book. 

Our results extend previous analysis of heavy traffic limits for such auction processes (Kruk 
(2003), Bayraktar et al. (2006), Cont and de Larrard (2010)) to a setting which is relevant and 
useful for quantitative modeling of limit order books and provide a foundation for recently proposed 
diffusion models for order book dynamics Avcllancda ct al. (2011). 

Outline. The paper is organized as follows. Section 2 describes a general framework for the 
dynamics of a limit order book; various examples of models studied in the literature are shown to 
fall within this modeling framework (Section 2.3). Section 3 reviews some statistical properties of 
high frequency order flow in limit order markets: these properties highlight the complex nature of 
the order flow and motivate the statistical assumptions used to derive the diffusion limit. Section 4 
contains our main result: Theorem 2 shows that, in a limit order market where orders arrive at high 
frequency, the bid and ask queues behaves like a Markov process in the positive quadrant which 
diffuses inside the quadrant and jumps to the interior each time it hits the boundary. We provide 
a complete description of this process, and use it to derive, in Section 4.3, a simple jump- diffusion 
approximation for the joint dynamics of bid and ask queues, which is easier to study and simulate 
than the initial queueing system. 

In particular, we show that in this asymptotic regime the price process is characterized as a 
piecewise constant process whose transition times correspond to hitting times of the axes by a two 
dimensional Brownian motion in the positive orthant (Proposition 1). This result allows to study 
analytically various quantities of interest, such as the distribution of the duration between price 
moves and the probability of an increase in the price: this is discussed in Section 5. 

2. A model for the dynamics of a limit order book 
2.1. Reduced-form representation of a limit order book 

Empirical studies of limit order markets suggest that the major component of the order flow occurs 
at the (best) bid and ask price levels (see e.g. Biais et al. (1995)). All electronic trading venues also 
allow to place limit orders pegged to the best available price (National Best Bid Offer, or NBBO); 
market makers used these pegged orders to liquidate their inventories. Furthermore, studies on the 
price impact of order book events show that the net effect of orders on the bid and ask queue 
sizes is the main factor driving price variations (Cont ct al. (2010a)). These observations, together 
with the fact that queue sizes at the best bid and ask of the consolidated order book are more 
easily obtainable (from records on trades and quotes) than information on deeper levels of the 
order book, motivate a reduced-form modeling approach in which we represent the state of the 
limit order book by 

• the bid price and the ask price s° 

• the size of the bid queue representing the outstanding limit buy orders at the bid, and 

• the size of the ask queue representing the outstanding limit sell orders at the ask 
Figure 3 summarizes this representation. 

If the stock is traded in several venues, the quantities and q"^ represent the best bids and offers 
in the consolidated order book, obtained by aggregating over all (visible) trading venues. At every 
time t, q\ (resp. q^) corresponds to all visible orders available at the bid price (resp. s") across 
all exchanges. 

The state of the order book is modified by order hook events: limit orders (at the bid or ask), 
market orders and cancelations (see Cont ct al. (2010b, a), Smith et al. (2003)). A limit buy (resp. 
sell) order of size x increases the size of the bid (resp. ask) queue by x, while a market buy (resp. 
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Figure 3 Simplified representation of a limit order book. 

sell) order decreases the corresponding queue size by x. Cancelation of x orders in a given queue 
reduces the queue size by x. Given that we are interested in the queue sizes at the best bid/ask 
levels, market orders and cancelations have the same effect on the queue sizes {q^,q^). 

The bid and ask prices are multiples of the tick size S. When either the bid or ask queue is 
depleted by market orders and cancelations, the price moves up or down to the next level of 
the order book. The price processes are thus piecewise constant processes whose transitions 

correspond to hitting times of the axes {(0, y), y > 0} U {(x, 0), x > 0} by the process qt = {q^, 

If the order book contains no 'gaps' (empty levels), these price increments are equal to one tick: 

• when the bid queue is depleted, the (bid) price decreases by one tick. 

• when the ask queue is depleted, the (ask) price increases by one tick. 

If there are gaps in the order book, this results in 'jumps' (i.e. variations of more than one tick) in 
the price dynamics. We will ignore this feature in what follows but it is not hard to generalize our 
results to include it. 

The quantity s° — is the bid-ask spread, which may be one or several ticks. As shown in Table 
3, for liquid stocks the bid-ask spread is equal to one tick for more than 98% of observations. 



Bid-ask spread 


1 tick 


2 tick 


> 3 tick 


Citigroup 


98.82 


1.18 





General Electric 


98.80 


1.18 


0.02 


General Motors 


98.71 


1.15 


0.14 



Table 3 Percentage of observations with a given bid-ask spread (June 26th, 2008). 

When either the bid or ask queue is depleted, the bid-ask spread widens immediately to more 
than one tick. Once the spread has increased, a flow of limit sell (resp. buy) orders quickly fills 
the gap and the spread reduces again to one tick. When a limit order is placed inside the spread, 
all the limit orders pegged to the NBBO price move in less than a millisecond to the price level 
corresponding to this new order. Once this happens, both the bid price and the ask price have 
increased (resp. decreased) by one tick. 

The histograms in Figure 4 show that this 'closing' of the spread takes place very quickly: as 
shown in Figure 4 (left) the lifetime of a spread larger than one tick is of the order of a couple of 
milliseconds, which is negligible compared to the lifetime of a spread equal to one tick (Figure 4 , 
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right). In our model we assume that the second step occurs infinitely fast: once the bid-ask spread 
widens, it returns immediately to one tick. For the example of Dow Jones stocks (Figure 4 ), this is 
a reasonable assumption since the widening of the spread lasts only a few milliseconds. This simply 
means that we are not trying to describe/model how the orders flow inside the bid- ask spread at 
the millisecond time scale and, when we describe the state of the order book after a price change 
we have in mind the state of the order book after the bid-ask spread has returned to one tick. 




1 2 3 4 5 B 0.5 1 1.6 2 2,5 3 

I lo' 

Figure 4 Left: Average lifetime, in milliseconds of a spread larger than one tick for Dow Jones stocks. Right: 
Average lifetime, in milliseconds of a spread equal to one tick. 

Under this assumption, each time one of the queues is depleted, both the bid queue and the 
ask queues move to a new position and the bid-ask spread remains equal to one tick after the 
price change. Thus, under our assumptions the bid-ask spread is equal to one tick, i.e. = 6, 
resulting in a further reduction of dimension in the model. 

Once either the bid or the ask queue are depleted, the bid and ask queues assume new values. 
Instead of keeping track of arrival, cancelation and execution of orders at all price levels (as in 
Cont et al. (2010b), Smith et al. (2003)), we treat the queue sizes after a price change as a 
stationary sequence of random variables whose distribution represents the depth of the order book 
in a statistical sense. More specifically, we model the size of the bid and ask queues after a price 
increase by a stationary sequence {Rk)k>i of random variables with values in N^. Similarly, the size 
of the bid and ask queues after a price decrease is modeled by a stationary sequence {Rk)k>i of 
random variables with values in N^. The sequences {Rk)k>i and {Rk)k>i summarize the interaction 
of the queues at the best bid/ask levels with the rest of the order book, viewed here as a 'reservoir' 
of limit orders. 

The variables R^ (resp. Rk) have a common distribution which represents the depth of the order 
book after a price increase (resp. decrease): Figure 5 shows the (joint) empirical distribution of bid 
and ask queue sizes after a price move for Citigroup stock on June 26th 2008. 

The simplest specification could be to take {Rk)k>i, {Rk)k>i to be IID sequences; this approach, 
used in Cont and de Larrard (2010), turns out to be good enough for many purposes. But this IID 
assumption is not necessary; in the next section we will see more general specifications which allow 
for serial dependence. 

In summary, state of the limit order book is thus described by a continuous-time process 
(S(,(7(,g") which takes values in the discrete state space (5Z x N^, with piecewise constant sample 
paths whose transitions correspond to the order book events. Denoting by (t°,i > 1) (resp. i^) the 
event times at the ask (resp. the bid), V"" (resp. V^'^) the corresponding change in ask (resp. bid) 
queue size, and k{t) the number of price changes in [0,t], the above assumptions translate into the 
following dynamics for {s^,q^,q1'): 
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Figure 5 Joint density of bid and ask queue sizes after a price move (Citigroup, June 26th 2008). 

• If an order or cancelation of size arrives on the ask side at t = t°;, 

— if + V^'^ > 0, the order can be satisfied without changing the price; 

— if + <0, the ask queue is depleted, the price increases by one 'tick' of size 5, and the 
queue sizes take new values Rk{t) = {Rk{t)^-^k{t))^ 

{sUlq-i) = {s\_,qUqt-+ynU_>-vt + (^L + 5, (1) 

• If an order or cancelation of size arrives on the bid side at t = t\^ 

— if q'l_ + V^** > 0, the order can be satisfied without changing the price; 

— if q\_ + < 0, the bid queue gets depleted, the price decreases by one 'tick' of size 5 and 
the queue sizes take new values Rk{t) = {^X(t)i^k(t))'- 

(slqlqt) = {s\_,qt+V,\qt_)\,_^_y, + {s\_ - 5,Rl^,^,Rl(,^)\,_^_y.. (2) 

2.2. The limit order book as a 'regulated' process in the orthant 

As in the case of reflected processes arising in queueing networks, the process qt = [ql^qt) niay be 
constructed from the net order flow process 

where (resp. N^) is the number of events (i.e. orders or cancelations) occurring at the bid (resp. 
the ask) during [0,t]. Xt = (xt,x°) is analogous to the 'net input' process in queuing systems Whitt 
(2002): x^ (resp. x^) represents the cumulative sum of all orders and cancelations at the bid (resp. 
the ask) between and t. 

Q — {l\i(lt)t>o which takes values in the positive orthant, may be constructed from x by reini- 
tializing its value to a a new position inside the positive orthant according to the rules (l)-(2) 
each time one of the queues is depleted: every time {qt)t>o attempts to exit the positive orthant, 
it jumps to a a new position inside the orthant, taken from the sequence (i?„,i?„). 

This construction may be done path by path, as follows: 
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Definition 1. Let w G Z)([0,oo),]R^) be a right-continuous function with left hmits (i.e. a cadlag 
function), R = (i?„)„>i and R = {Rn)n>i two sequences with values in R^. There exists a unique 
cadlag function '^{u,R,R) G D([0,oo),M^) such that 

• For t<ri, let ^{uj,R,R){t) = uj{t) where 

Ti = inf{t > 0, w(t).(l,0) < or a;(t).(0, 1) < 0}. 

is the first exit time of cj from the positive orthant. 

• ^{u;,R,R){ti) = Ri if ^{uj, R, R){ti-).{0,1) < 0, and ^{oj, R, R){ti) = Ri if 
^{uj,R,R){n-).{l,0)<0. 

• For fc > 1: 

^{uj, R, R){t + Tk) = ^-(0;, R, R){n) + uj{t + Tk) - uj{Tk) for < f < r^+i - r^, where 

Tk = \ni{t>Tu-i, a;(t).(l,0)<^(L^,i?,i?)(rfc_i)-(l,0) oi u{t).{Q,l) <-^{u,R,R){Tk-i).{^A)] 

• ^{u,R,R){n) = Rkif^{u;,R,R){Tk-).{0^)<^{uj,R,R){Tk-i).{l,0)and^{uj,R,R) 
otherwise. 

The path ^'(a;, R, R) is obtained by "regulating" the path uj with the sequences {R, R): in between 
two exit times, the increments of ^{uj, R, R) follow those of uj and each time the process attempts 
to exit the positive orthant by crossing the x-axis (resp. the y-axis), it jumps to a a new position 
inside the orthant, taken from the sequence (i?„)„>i (resp. from the sequence (-R„)„>i). 

Unlike the more familiar case of a continuous reflection at the boundary, which arises in heavy- 
traffic limits of multiclass queueing systems (see Harrison (1978), Harrison and Nguyen (1993), 
Whitt (2002), Ramanan and Rciman (2003) for examples), this construction introduces a disconti- 
nuity by pushing the process into the interior of the positive orthant each time it attempts to exit 
from the axes. 

To study the continuity properties of this map, we endow Z?([0, oo),M^) with Skorokhod's Ji 
topology Billingsley (1968), Lindvall (1973) and the set (M^) with the topology induced by 'cylin- 
drical' semi-norms, defined as follows: for a sequence (i?")„>i in (M^) 

i?"™i?G (M^)" ^ (yk>l, sup{|i??^-i?i|,...,|i?l!-i?fe|)™0) . 

D([0,oo),M2) X [Rlf X (Rl)" is then endowed with the corresponding product topology. 

Theorem 1. Let R = (i?„)„>i,i?= (i?„)„>i be sequences in ]0, oo[x]0, oo[ which do not have any 
accumulation point on the axes. If to £ C°([0, 00), M^) is such that 

iO,0)^^{u;,R,R)i[0,oo)). (3) 

Then the map 

'^■.D{[0,oo),R'') X {Rlf X {Rlf D([0,oo),M^) (4) 

is continuous at (uj,R,R). 

Proof: see Section 6.2 in the Appendix. 

This construction may be applied to any cadlag stochastic process: given a cadlag process X 
with values in and (random) sequences R = {R,i)n>i and R = {Rn)n>i with values in R"^, the 
process \I'(X, i?, i?) is a cadlag process with values in M^. 

It is easy to see that the order book process qt = {q^, q1) may be constructed by this procedure: 

Lemma 1. q = {q^ ^q"-) = ^{x,R,R) where 
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• Xt = {x\,x1) = (^XliJi ^i^^ili ^tj order flow at the hid and the ask, 

• R = {Rn)n>i is the sequence of queue sizes after a price increase, and 

• R = {Rn)n>i is the sequence of queue sizes after a price decrease. 

One can thus build a statistical model for the limit order book by specifying the joint law of x 
and of the regulating sequences {R,R). This approach simplifies the study of the (asymptotic) 
properties of qt = (q^jQl). 

Example 1 (IID reinitializations). The simplest case is the case where the queue length after 
each price change is independent from the history of the order book, as in C'oiit and dc Larrard 
(2010). R = (i?„)„>i and R = (-R„)„>i are then IID sequences with values in ]0, oo[^. Figure 5 shows 
an example of such a distribution for a liquid stock. 

The law of the process Q = ^{x, R, R) is then entirely determined by the law of the net order 
flow x and the distributions of -R„, it can be constructed from the concatenation of the laws 
of (xf , Tfe <t< Tfc+i) for A; > (where we define Tq := 0). 

Example 2 (Pegged limit orders). Most electronic trading platforms allow to place hmit or- 
ders which are pegged to the best quote: if the best quote moves to a new price level, a pegged limit 
order moves along with it to the new price level. The presence of pegged orders leads to positive 
autocorrelation and dependence in the queue size before/after a price change. The queue size after 
a price change may be modeled as 

• qr„ = Rn = {e^ + Pl^^-^i-n) ^he price has increased, and 

• qr^ = Rn = i^n^^n + Pirn") pricB has dccreascd 

where e„ = (e(^,e°),e„ = (e^,e°) are IID sequences. Empirically, one observes a correlation of ~ 
10% — 20% between the queue lengths before and after a price change, which suggests an order 
magnitude for the fraction of pegged orders. 

As in the previous example, the law of of the process q = ^{x, R, R) is determined by the law of 
the net order flow x, the coefficients /3,/3 and the distributions of e,e: it can be constructed from 
the concatenation of the laws of {xt,Tk <t< Tk+i) for A; > 0. 

More generally, one could consider other extensions where the queue size after a price move may 
depend in a (nonlinear) way on the queue size before the price move and a random term e„ 
representing the inflow of new orders after the n-th price change: 

9r„ =5(9r„-,e„). (5) 

The results given below hold for this general specification although the examples 1 and 2 above 
are sufficiently general for most applications. 

2.3. Examples 

The framework described in Section 2.1 is quite general: it allows a wide class of specifications 
for the order fiow process, and contains as special cases various models proposed in the literature. 
Each model involves a specification for the (random) sequences {f^ ,t^,V^'^ ,V^)i>i, R = {Rn)n>i and 
R = {Rn)n>i or, equivalents, (I;^ V^,V^),>i, R = (i?„)„>i and R = (i?„)„>i where T^^ = t^^, - 1^ 
(resp. T^' = — t^) are the durations between order book events on the ask (resp. the bid) side. 

2.3.1. Models based on Poisson point processes Cont and dc Larrard (2010) study a 
stylized model of a limit order market in which market orders, limit orders and cancelations arrive 
at independent and exponential times with corresponding rates ^, A and 9, the process q = {q^, g°) 
becomes a Markov process. If we assume additionally that all orders have the same size, the 
dynamics of the reduced limit order book is described by: 
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• The sequence (T°)i>o is a sequence of independent random variables with exponential distri- 
bution with parameter X + 6 + /j,, 

• The sequence {Tl')i>Q is a sequence of independent random variables with exponential distri- 
bution with parameter X + 6 + ^, 

• The sequence {V^)i>o is a sequence of independent random variables with 

nvr = 1] = ^nd nvt = -1] ^ + ^ 



The sequence {V^)i>o is a sequence of independent random variables with 

= 1] = , / , , and ¥[¥!> = -1] 



• All these sequences are independent. 
It is readily verified that this model is a special case of the framework of Section 2.1: (qt) t>o may 
be constructed as in Definition 1, where the unconstrained process Xt is now a compound Poisson 
process. 

2.3.2. Self-exciting point processes Empirical studies of order durations highlight the de- 
pendence in the sequence of order durations. This feature, which is not captured in models based 
on Poisson processes, may be adequately represented by a multidimensional self-exciting point pro- 
cess Andersen et al. (2010), Hautsch (2004), in which the arrival rate Xi{t) of an order of type i 
is represented as a stochastic process whose value depends on the recent history of the order flow: 
each new order increases the rate of arrival for subsequent orders of the same type (self-exciting 
property) and may also affect the rate of arrival of other order types (mutually exciting property): 



A,(i) = + ^S,, [ e-'^^^'-'UN,{s) 
Jo 



i=i 

Here 5ij measures the impact of events of type j on the rate of arrival of subsequent events of type 
i: as each event of type j occurs, A; increases by 6ij. In between events, Xi{t) decays exponentially 
at rate k^. Maximum likelihood estimation of this model on TAQ data Andersen cl al. (2010) shows 
evidence of self-exciting and mutually exciting features in order flow: the coefficients (5^ are all 
significantly different from zero and positive, with > Sij for j ^ i. 

2.3.3. Autoregressive conditional durations Models based on Poisson process fail to cap- 
ture serial dependence in the sequence of durations, which manifests itself in the form of clustering 
of order book events. One approach for incorporating serial dependence in event durations is to 
represent the duration between transactions i — 1 and i as 

where (ei)i>i is a sequence of independent positive random variables with common distribution 
and E[ei] = 1 and the conditional duration tpi = K[Ti\'ipi^j,Ti^j,j > 1] is modeled as a function of 
past history of the process: 

= G('i/'j-l,'0*-2, .■■,T,^i,T,_2, ••■,••)■ 

Engle and Russell's Autoregressive Conditional Duration model Engie and Russell (1998) propose 
an ARMA(p, q) representation for G: 



-0, = flo + ^ aktpi-k + ^ bgTj^k 



Cont, Larrard: Order book dynamics in liquid markets 



13 



where (ao,...,ap) and {bi,...,bg) are positive constants. The ACD-GARCH model Ghysels and 
Jasiak (1998) combine this model with a GARCH model for the returns. Engie (2000) proposes 
a GARCH-type model with random durations where the volatility of a price change may depend 
on the previous durations. Variants and extensions are discussed in Hautscli (2004). Such models, 
like ARMA or GARCH models defined on fixed time intervals, have likelihood functions which are 
numerically computable. Although these references focus on transaction data, the framework can 
be adapted to model the durations (T°,i > 1) and {Tj^ji > 1) between order book events with the 
ACD framework Hautsch (2004). 

2.3.4. A limit order market with patient and impatient agents Another way of spec- 
ifying a stochastic model for the order flow in a limit order market is to use an 'agent-based' 
formulation where agent types are characterized in terms of the statistical properties of the order 
flow they generate. Consider for example a market with three types of traders: 

• impatient traders who only submit market orders: 

• patient traders who use only limit orders: this is the case for example of traders who place 
stop loss orders or engage in strategies such as mean-reversion arbitrage or pairs trading which are 
only profitable with limit orders. 

• other traders who use both limit and market orders; we will assume these traders submit a 
proportion 7 of their orders as limit orders and (1 — 7) as market orders, where < 7 < 1. 
Denote by m (resp. I) the proportion of orders generated by impatient (resp. patient) traders: 

Vi > 1, P[i — th trader uses only market orders] = m, 
P[i — th trader uses only limit orders] = I, 
P[ith trader uses both limit and market orders] = 1 — I — m. 

Assume that the sequence {Ti,i > 1) of duration between consecutive orders is a stationary ergodic 
sequence of random variables with E[Ti] < 00, that each trader has an equal chance of being a 
buyer or a seller and that the type of trader (buyer or seller) is independent from the past: 

P[i — th trader is a buyer] = P[i — th trader is a seller] = - 

Trader i generates an order of size V^, where (T^,i > 1) is an IID sequence with: 

777 

mv!', vn = {v,M = p[(y/, vn = (o, v,)] = -, 
vn = i-y^M = p[(y/, vn = (o, -v,)] = ^, 

1 — / — 777 

nivi',vn = m,-a-im]=n{v^vn = {-{^-im,iy)] = — - — . 

3. Statistical properties of high-frequency order flow 

As described in Section 2.1, the sequence of order book events -the order flow- is characterized by 
the sequences (T", i > 1) and {T^, z > 1) of durations between orders and the sequences of order sizes 
> 1) and {V^,i > 1). In this section we illustrate the statistical properties of these sequences 
using high-frequency quotes and trades for liquid US stocks -CitiGroup, General Electric, General 
Motors- on June 26th, 2008. 
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3.1. Order sizes 

Empirical studies Bouchaud et al. (2002, 2008), Gopikrishnan et al. (2000), Maslov and Mills (2001) 
have shown that order sizes are highly heterogeneous and exhibit heavy-tailed distributions, with 
Pareto-type tails: 

P{V," > x) ~ Cx"^ 

with tail exponent /3 > between 2 and 3, which corresponds to a series with finite variance but 
infinite moments of order > 3. The tail exponent /3 > is difficult to estimate precisely, but the 
Hill estimator Resnick (2006) can be used to measure the heaviness of the tails. Table 4 gives the 
Hill estimator of the tail coefficient of order sizes for our samples. This estimator is larger than 2 
for both the bid and the ask; this means that the sequence of order sizes have a finite moment of 
order two. 





Bid side 


Ask side 


Citigroup 


[0.42, 0.46] 


[0.29, 0.32] 


General Electric 


[0.42, 0.45] 


[0.41, 0.46] 


General Motors 


[0.36, 0.42] 


[0.44, 0.51] 



Table 4 95-percent confidence interval of the Hill estimator of the sequence of order sizes. When the Hill 
estimator is < 0.5, the estimated tail index is large than 2 and the distribution has finite variance. 



The sequences of order sizes (V"j°,i > 1) and {V^,i > 1) exhibit insignificant autocorrelation, as 
observed on Figure 6. However, they are far from being independent: the series of squared order 
sizes {{V^y,i > 1) and {{V^y,i > 1) are positively correlated, as revealed by their autocorrelation 
functions (displayed in Figure 7). 



Sample Autocorrelation Function (ACF) 



Sample Autocorrelation Function (ACF) 
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Figure 6 Autocorrelogram of the sequence of order sizes. Order coming at the ask on the left and at the bid on 
the right. 



Finally, the sequences {V°',i > 1) and (T^^,i > 1) may be negatively correlated. This stems from 
the fact that a buyer can simultaneously use market orders on the ask side (which correspond to 
negative values of 1^° and limit orders on the bid side (which correspond to positive values of V^^); 
the same argument holds for sellers (see Section 2.3.4). 
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Sample Autocorrelation Function (ACF) 
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Sample Autocorrelation Function (ACF) 
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Figure 7 Autocorrelogram of the sequence of absolute order sizes. Order coming at the ask on the left and at the 
bid on the right. 



These properties of the sequence (V^ ,V^)i>i may be modeled using a bivariate ARCH process: 

{a':)' = a'o + al {V^f, {a^f = al + J^ where {z^ zt)..>, N (^0, (^J ) 

and (oq, a^, Oq, a") are positive coefficients satisfying 

0<al + al<l, and 0<< + <<l. (6) 

As shown by Bougcrol and Picard (1992), under the assumption (6), the sequences of order sizes 
(V"/, ^ > 1) and {V°',i > 1) is then a weh defined, stationary sequence of random variables with finite 
second-order moments, satisfying the properties enumerated above. 

3.2. Durations 

The timing of order book events is describe by the sequence of durations {Tj',i > 1) at the bid 
and iT°',i > 1) at the ask. These sequences have zero autocorrelation (see Figure 8) but are not 
independence sequences: for example, as shown in Figure 9, the sequence of inverse durations 
(l/Tf,? > 1) and {1/Tl',i > 1) is strongly correlated in each case. 



Figure 10 represents the empirical distribution functions P[T° > u]) and F[T^ > u] in logarithmic 
scale. Both empirical distributions exhibit thin, exponential-type tails (which implies in particular 
that and have finite expectation). 
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Sample Autocorrelation Function (ACF) 



Sample Autocorrelation Function (ACF) 
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Figure 8 Autocorrelogram of the sequence of durations for events at the ask (left) and the bid (right). 



Sample Autocorrelation Function (ACF) 
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Figure 9 Autocorrelogram of the sequence of inverse durations for events at the ask (left) and the bid (right). 
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4. Heavy Traffic limit 

At very high frequency, the limit order book is described by a two-dimensional piecewise constant 
process qt = (g^ , (?")t>o, whose evolution is determined by the flow of orders. The complex nature 
of this order flow -heterogeneity and serial dependence in order sizes, dependence between orders 
coming at the ask and at the bid- described in section 3, makes it difficult to describe qt in an 
analytically tractable manner which would allow the quantities of interest to be computed either in 
closed form or numerically in real time applications. However, if one is interested in the evolution 
of the order book over time scales much larger than the interval between individual order book 
events, the (coarse-grained) dynamics of the queue sizes may be described in terms of a simpler 
process Q, called the heavy traffic approximation of q. In this limit, the complex dynamics of the 
discrete queueing system is approximated by a simpler system with a continuous state space, which 
can be either described by a system of ordinary differential equations (in the 'fluid limit', where 
random fluctuations in queue sizes vanish) or a system of stochastic differential equations (in the 
'diffusion limit' where random fluctuations dominate). This idea has been widely used in queueing 
theory to obtain useful analytical insights into the dynamics of queueing systems Harrison and 
Nguyen (1993), Iglchart and Whitt (1970), Whitt (2002). 

We argue in this section that the heavy traffic limit is highly relevant for the study of limit order 
books in liquid markets, and that the correct scaling limit for the liquid stocks examined in our 
data sets is the "diffusion" limit. This heavy traffic limit is then derived in Section 4 and described 
in Section sec.Markovianapproximation. 

4.1. Fluid limit or diffusion limit? 

One way of viewing the heavy traffic limit is to view the limit order book at a lower time resolution, 
by grouping together events in batches of size n. Since the inter-event durations are finite, this is 
equivalent to rescaling time by n. The impact, on the net order flow, of a batch of n events at the 
ask is 

v{^+v,^+v,^ + ...+v: _ iv,^-W) + iv,^-W) + ... + iv:-W) 

I y ^ y 5 



n yjn 



where (y^ > 1) is the sequence of order sizes at the ask and V'^ = E[V"]. Under appropriate 
assumptions (see next section), this sum behaves approximately as a Gaussian random variable for 
large n: 

"1 ^''2 ^"3 ^•••.^Ir^_^(^ya^Var(yi")) as n^oo. (7) 



n 



Two regimes are possible, depending on the behavior of the ratio J^yg^ as n grows: 

• If /" — )• oo as n — )• oo, the correct approximation is given by the fluid limit, which describe 

^Var(lA") 

the (deterministic) behavior of the average queue size. 

• If lim„^oo /" < oo, the rescaled queue sizes behave like a diffusion process. 

The fluid limit corresponds to the regime of law of large numbers, where random fluctuations 
average out and the limit is described by average queue size, whereas the diffusion limit corre- 
sponds to the regime of the (functional) central limit theorem, where fluctuations in queue size are 

asymptotically Gaussian. 

Figure 11 displays the histogram of the ratio /" for stocks in the Dow Jones index, where for 

each stock n is chosen to represent the average number of order book events in a 10 second interval 
(typically n ~ 100 — 1000). This ratio is shown to be rather small at such intraday time scales, 
showing that the diffusion approximation, rather than the fluid limit, is the relevant approximation 
to use here. 
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Figure 11 Empirical distribution of the ratio y/n/V" ^Var(V]") showing the relative importance of average change 
vs fluctuations in queue size, for stocks in the Dow- Jones index during June 08 (see Section 4.1). Low 
values of the ratio indicate that intraday changes in bid/ask queue size are dominated by fluctuations, 
rather than the average motion of the queue. Left: bid side. Right: ask side. 

Indeed, bid and ask queue sizes [q^^q^) exhibit a diffusion-type behavior at such intraday time 
scales: Figure 12 shows the path of the net order flow process 

/ A^t" JVf \ 

Xt = {qU^,)+\Y.V^,Y.^t\ (8) 
\ i=i i=i / 

sampled every second for CitiGroup stocks on a typical trading day. In this example, for which the 
average time between consecutive orders is A^^ ~ 13 ms ^ 1 second, we observe that the process 
X behaves like a diffusion in the orthant with negative drift: the randomness of queue sizes does 
not average out at this time scale. 




Figure 12 Evolution of the net order flow Xt — {X^^Xf) given in Eq. (8) for CitiGroup shares over one trading 
day (June 26, 2008). 
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We will now show that this is a general result: under mild assumptions on the order flow process, 
we will show that the (rescaled) queue size process 

(^,^),>o (9) 



converges in distribution to a Markov process (Qt)t>o in the positive orthant, whose features we 
will describe in terms of the statistical properties of the order flow. 

4.2. A functional central limit theorem for the limit order book 

Consider now a sequence = ((/")t>o of processes, where represents the dynamics of the bid and 
ask queues in the limit order book at a time resolution corresponding to n events (see discussion 
above). The dynamics of is characterized by the sequence of order sizes (T^"''', VJ"''^)i>i, durations 
(T"'^, rj"''")i>i between orders and the fact that, at each price change 

• g"^ = R^ = g{q^^_,e^^) if the price has increased, and 

• g"^ = i?^ = eJJ) if the price has decreased, 

where {€^,k > 1) is an IID sequence with distribution /„, and {€]^,k > 1) is an IID sequence with 
distribution /„. Note that this specification includes Examples 1 and 2 as special cases. 

We make the following assumptions, which allow for an analytical study of the heavy traffic limit 
and are sufficiently general to accommodate high frequency data sets of trades and quotes such as 
the ones described in Section 3: 

Assumption 1. {r^-'^ ^t'^'") i>i is a stationary array of positive random variables whose common 
distribution has a continuous density and satisfies 

hm — — = < oo, hm — — = < oo. 

ra-i-oo n A" n->oo n A 

A" (resp. A'') represents the arrival rate of orders at the ask (resp. the bid). 

Assumption 2. V^"'^)i>i is a stationary, uniformly mixing array of random variables satis- 

fying 

V^E[yi"''^]™F^, V^E[Vi"^'']™y^, (10) 

oo 

lim E[(y,"''^ - F^)2] + 2 V coY{Vr \ l^"'") = vl<oo, and 

oo 

lim ¥.[{V^^' - V~bf] + 2 V cov(yl"'^ V^^') = vl < oo. 

1=2 

The assumption of uniform mixing (Billingsley 1968, Ch. 4) implies that the partial sums of order 
sizes verify a central limit theorem, but allows for various types of serial dependence in order sizes. 
The scaling assumptions on the first two moments corresponds to the properties of the empirical 
data discussed in Section 4.1. Under Assumption 2, one can define 

p:= lim — ( 2max(A^A')cov(yl"'^yl"■') + 2y AW(yl"^^y,"'') + A%v(yl"'^y,"■") I . (11) 



1=1 



p£ (—1,1) may be interpreted as a measure of 'correlation' between event sizes at the bid and 
event sizes at the ask. 

These assumptions hold for the examples of Section 2.3. In the case of the Hawkes model. 
Assumption 1 was shown to hold in Bacry et al. (2010). Also, these assumptions are quite plausible 
for high frequnecy quotes for liquid US stocks since, as argued in Section 3: 
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• The tail index of order sizes is larger than two, so the sequences {V^^,i > 1) and {V^^i > 1) 
have a finite second moment. 

• The sequence of order sizes is uncorrelated i.e. has statistically insignificant autocorrelation. 
Therefore the sum of autocorrelations of order sizes is finite (zero, in fact). 

• The sequence of inter-event durations has a finite empirical mean and is not autocorrelated. 
These empirical observations support the plausibility of Assumptions 1 and 2 for the data sets 
examined. 

Assumption 2 has an intuitive interpretation: if orders are grouped in batches of n orders, then 
Assumption 2 amounts to stating that the variance of batch sizes should scale linear with n. This 
assumption can be checked empirically, using a variance ratio test for example: Figure 13 shows 
that this linear relation is indeed verifies for the data sets examined in Section 3. 




Figure 13 Variance of batch sizes of n orders, for General Electric shares, on June 26th, 2008. Left: ask side. 
Right: bid side. 



The following scaling assumption states that, when grouping orders in batches of n orders, a good 
proportion of batches should have a size 0{y/n) (otherwise their impact will vanish in the limit 
when n becomes large): 

Assumption 3. There exist probability distributions F,F on the interior (0,oo) x (0,oo) of the 
positive orthant, such that 



nfniVn .) ^ F and n/„(A/n .) =^> F. 
Assumption 4. ^ g C^(M^ x ]R^,]0, oop) and 

3a>0,V(x,y) gR^ xR^, \\g{x,y)\\> amm{yi,y2). 
Finally, we add the following condition for the initial value of the queue sizes: 

/ n.b n,a \ 

iQo qo -.-^(^^,y^)-,(o,o). (12) 



n vn 



The following theorem, whose proof is given in the Appendix, describes the joint dynamics of the 
bid and ask queues in this heavy traffic limit: 
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Theorem 2 (Heavy traffic limit). Under Assumptions 3, 1, 2 and 4, the rescaled process 

{Q^,t>0)=(%,t>0 



n 

2 



converges weakly, on the Skorokhod space (D(]0,oo[,]R^), , Ji), 

to a Markov process {Qt)t>o with values in — {(0,0)} and infinitesimal generator Q given by 
^, / N , ^r^dh , TTzdh Xavl d'^h Xtv? d'^h /- — — d'^h , . 

eft(i,0)=/ ghig{{x,0),(u,v)))F{du,dv), gh(0,y)=[ gh(gmy),(u,v)))F{du.dv), (U) 

for X > 0,y > 0, whose domain is the set dom{Q) of functions /i G C^(]0, oo[x]0, oo[,]R) n C"(]R^,]R) 
verifying the boundary conditions 



Vx>0, h{x,0)= / h{g{{x,0),{u,v)))F{du,dv) = 0, 
yy>0, h{0,y)= [ h{g{{0,y),{u,v)))F{du,dv) = 0. 



Proof. We outline here the main steps of the proof. The technical details are given in the 
Appendix. Define the counting processes 

Nr = sup{fc > 0, Tf^" + ... + T,"'" < t} and iVf'" = sup{fc > 0, + ... + T,'^" < t} (15) 

which correspond to the number of events at the ask (resp. the bid), and the net order flow 

/ nt -iTh,n 7it Tra.n 



Jn ^ — ' \/n 

i=l ^ 1=1 



Then, as shown in Proposition 3 (see Appendix), X" converges in distribution on (Z?([0, cxd[,M^), Ji] 
to a two-dimensional Brownian motion with drift 



(xnt>o ^ ^^t+^(A''y^A"y-; 

where Z is a planar Brownian motion with covariance matrix 



t>0 



Vb 



Under assumption 3, using the Skorokhod representation theorem, there exist IID sequences 
((eJJ, n > 1), (e^, n > l),ek,ik)k>i and a copy X of the process 



{xo,yo) + Zt + t{X''V\X''V-: 



t>o 



on some probability space {^lo,13, Q) such that e^' ~ /„, ^ fn,ek ^ F,ek ^ F and 



X" "=r X ; VA; > 1, ^ "4°° e,, ^ ™ ik ] = 1. 



Using the notations of Appendix 6.2, denote by 
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• r" =r(X") the first exit time of X" from the interior ]0, oo[x]0, oo[ of the orthant and 

• tJ! the first exit time of ^'fe(X", Q;^„, Q;^n_ ) from ]0, oo[x]0, oo[. 

We can now construct the process Q by an induction procedure. Let Ti = t{X) be the first exit time 
of X from the orthant. By definition of Q, Qt = Xt for t <Ti and, by continuity of the first-passage 
time map and the last-evaluation map at a first passage time (Whitt 2002, Sec. 13.6.3), 

(rr,Q:!„_)™(ri,Q.,_)Q-a.s. 

By definition of Q, we have 

Qti = g{Xri-,ei)lxr-^^.{o,i)<o + g{Xr-^_,ek)lx^-J^.{l,o)<o■ 
X is a continuous process and the probability thats its path crosses the origin is zero, so by Lemma 
2, X lies with probability 1 in the continuity set of the map G : cj — )• lc,j^(^j.(o,i)<o- So using the 
continuity of (/(., .), we can apply the continuous mapping theorem (Billingsley 1968, Theorem 5.1), 
to conclude that 

Q;;n"^°°Q,, q-a.s. 

Let us now show the induction step: assume that we have defined Q on [0,rfc_i] and shown that 

(ti", . . , rfc"_i , Q^^n , . . . , Q'^n_^) "4°°(ti ,..,n-uQr„:., Qr,,_, ) Q-a.s. 

Since Q((0,0) ^ ^k{X,Qr^, ...,Qr^_^){[0,oo)) ) = 1, Lemma 4 implies that {X,Qr^, ...,Qr^_J lies 
with probability 1 in the continuity set of ^'fc, so by the continuous mapping theorem 

^k{X\Q^n,...,Q",j.J''^°°^k{X,Qr„...,Qr,^,) Q-a.s. 

Define now as the first exit time of ^k{X,Qr^, ...,Q'^^__^) from ]0, oo[x]0, oo[. As before, by 
continuity of the first-passage time map and the last-value map at a first passage time (Whitt 
2002, Sec. 13.6.3), 

(r;,Q:!„_)™(rfc,Q.,_) Q-a.s. 

k 

We can now extend the definition of Q to [0,r/j] by setting 

Qt = '^k{X,Qr^,...,QrJ{t) for t<Tk, and 

Qrk =5(Qrfc-,efc)l*fc(X,Q^j,...,Q^j^_j).(0,l)<0 +fl'(<3rfc-,efc)l*fc(X,Q^-^,...,Q^j^_^).(l,0)<0 

As above, using the continuity properties of from Lemma 4 we conclude that (5"n — )• Qi-j. a-S. 
So finally, we have shown that 

VA:>1, (rr,..,r,",Q;!„,...,Q;;„)"4°°(Ti,..,T,,Q.,,...,g.J Q-a.s. 
We can now construct the sequences R, R by setting 

. Rk=Qr, ii'^k{X,Qr„...,Qr,_,){Tk-).{0,l)<0, 

. Rk = Qr,ii^k{X,Qr„...,Qr,_,){n-).{l,0)<0. 
Then Q = "^{X,R,R) where ^ is the map defined in Definition 1. Let us now show that {X,R,R) 
lies with probability 1 in the Ji— continuity set of ^, in order to apply the continuous mapping 
theorem. X is a continuous process whose paths lie in C"'([0, oo),]R^ — {(0, 0)}) almost surely. Since 
F and F have zero mass on the axes, with probability 1 the sequences (efc)fc>i, {ek)k>i do not have 
any accumulation point on the axes. Assumption 4 then implies that the sequences {Rk)k>i, iRk)k>i 
do not have any accumulation point on the axes. From the definition of ^ (Definition 1), Q jumps 



Cont, Larrard: Order book dynamics in liquid markets 



23 



at each hitting time of the axes and, in between two jumps, its increments follow those of the 
planar Brownian motion X. Since F, F have no mass at the origin and planar Brownian paths have 
a zero probability of hitting isolated points, with probability 1 the graph of Q = ^{X,R,R) does 
not hit the origin : 

q( (0,0)^M'(X,i?,i?)([0,oo) ) ) =1. (16) 

So the triplet {X, R, R) satisfies the conditions of Theorem 1 almost-surely i.e. ^ is continuous at 
{X,R,R) with probability 1. We can therefore apply the continuous mapping theorem (Billingsley 
1968, Theorem 5.1) and conclude that 

Q" = (X" , i?" , fl" ) Q = ^' (X, «) . 

The process Q = ^{X,R,R) can be explicitly construction from the planar Brownian motion X 
and the sequences R, R: Q follows the increments of X and is reinitialized to i?„ or i?„ at each 
hitting time of the axes. Lemma 5 in Appendix 6.4 uses this description to show that Q is a Markov 
process whose generator is given by (13)- (14). 

Remark 1 (Levy process limits). The diffusion approximation inside the orthant fails when 
order sizes do not have a finite second moment. For example, if the sequence {V°',V^'') is regularly 
varying with tail exponent a £ (0,2) (see Resnick (2006) for definitions), the heavy-traffic approx- 
imation Q is a pure-jump process in the positive orthant, constructed by applying the map ^' to a 
two-dimensional a-stable Levy process L: 

Q = -^{L,R,R), 

i.e. by re-initializing it according to (5) at each attempted exit from the positive orthant. We do 
not further develop this case here, but it may be of interest for the study of illiquid limit order 
markets, or those where order flow is dominated by large investors. 



4.3. Jump-difrusion approximation for order book dynamics 

Theorem 2 implies that, when examined over time scales much larger than the interval between 
order book events, the queue sizes q'' and are well described by a Markovian jump- diffusion 
process (Qt)t>o in the positive orthant ]R'j_ which behaves like a a planar Brownian motion with 
drift vector 

(A'y^, A^F^) (17) 

and covariance matrix 

in the interior ]0,oo[^ of the orthant and, at each hitting time of the axes, jumps to a new 
position 

• Q^^ = Rj^ = g{Qr^^,ek) if it hits the horizontal axis, 

• Q^^ = R^ = g[Q^^_,ek) if it hits the vertical axis, 

where the are IID with distribution F and the are IID with distribution F. We note that similar 
processes in the orthant were studied by Baccelli and Fayollc (1987) with queueing applications in 
mind, but not in the context of heavy traffic limits. 

This process is analytically and computationally tractable and allows various quantities related 
to intraday price behavior to be computed (see next section). 
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If 7o = (IE[Tf] +E[Tj''])/2 is the average time between order book events, (70 < 100 milliseconds), 
and 7i ^ 70 (typically, 71 ~ 10-100 seconds) then Theorem 2 leads to an approximation for the 
distributional properties of the queue dynamics in terms of Qt- 

qt^'^VN Qt/N where iV = ^ 

7o 

So, under Assumptions 3, 1, 2 and 4 the order book process {qt,Qt)t>a at the time scale 71 can be 
approximated by a Markov process which 

• behaves like a two-dimensional Brownian motion with drift (//;,, ^q) and covariance matrix A 
on {x>0}n{y> 0} with 

and, at at each hitting time of the axes, 

• jumps to a new value g{qt-, VNck) if q't^ = 0, 

• jumps to a new value g{qt^, VNck) if qt-^ = 0, 
where F, F are IID. 

This gives a rigorous justification for modeling the queue sizes by a diffusion process at such 
intraday time scales, as proposed in Avollancda ot al. (2U11). The parameters involved in this 
approximation are straightforward to estimate from empirical data: they involve estimating first 
and second moments of durations and order sizes. 

Example 3. Set for instance 71 = 30 seconds and 70 = {ElT"] +E[T''])/2. The following table 
shows the parameters (19) estimated from high frequency records or order book events for three 
liquid US stocks. 





Std deviation of 


Std deviation of 










Bid queue 


Ask queue 


Pb 




P 


Citigroup 


6256 


4457 


-1033 


-2467 


0.07 


General Electric 


2156 


2928 


-334 


-1291 


0.03 


General Motors 


578 


399 


-f78 


-96 


- 0.04 



Table 5 Parameters for the heavy-traffic approximation of bid / ask queues over a 30-second time scale. The 

unit is a number of orders per period of 30 seconds. 

In particular we observe that the order of magnitude of the standard deviation of queue lengths 
is an order of magnitude larger than their expected change. 

Example 4. Theorem 2 may also be used to derive jump- diffusion approximations for the limit 
order book in theoretical models such as the ones presented in Section 2.3. Let us illustrate this in 
the case of the heterogeneous trader model of Section 2.3.4. 

Let {T.i,i> 1) the sequence of duration between consecutive orders. We assume that this sequence 
is a sequence of stationary random variables with E[Ti] < 00. We also assume that every trader 
has an equal chance of being a buyer or a seller and that the type of trader (buyer or seller) is 
independent from the past: 

P[i — th trader is a buyer] = P[i — th trader is a seller] = - 

Finally the sequence of number of orders {Vi,i>l) is a stationary sequence of orders traded by 
the i-th trader with the property that IE[^^i'^] < 00. 
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This order flow given by {T^,i > 1), {Vi,i > 1), and the sequence of type (buyers or sellers, using 
limit orders, market orders or both) generates a sequence of durations {T°',i > 1), (T^;i > 1) and 
order sizes (V^",i > 1) and {Vf^,i > 1) which satisfy assumptions 1 and 2. 

The sequence of durations {T^,i > 1) and {Tj',i > 1) are two stationary sequences of random 
variables with finite mean: 



Vi>0, T, = T^ = T^. therefore E[r,] = E[i;"] = E[7;''] < oo. 
The sequence of order sizes {{V^ ^V""),! > 1) is a sequences of IID random variables with 

Vn = {V^M = Vt) = (0, m = -, (20) 

mv', v^n = i-v^M = mv', vn = (o, -v,)] = ^, (21) 

1 — / — m 

mv-, vn = m, - (i - m)] = mvn vn = (-(i - i)v,iv,)] = — - — . (22) 

Theorem 2 then shows that {Q^,Q"^) is a Markov process which behaves like a two-dimensional 
Brownian motion with drift (/ib,^^) and covariance matrix A inside the positive orthant {x > 
0} n {y > 0} where: 

/^6 = /ia = ^^(2m + 27(l-/-m)-l), A = v^(^y^, where (23) 

v^= ^ ' ^ m-t- +- ^ —(l-l-m)] and p = ) " , , <0. 

4 V 2 ^ J l + (l-/-m)(72-7-l/2) 

(24) 

Figure 14 displays the value of the correlation p in different scenarios as a function of 7 and the 
proportion 1 — (/ + m) of traders submitting orders of both types. 




Figure 14 Correlation p between bid and ask queue sizes for different scenario. 1 - (l + m) represents the proportion 
of traders using both market and hmit orders, 7 the proportion of limit orders and (1 — 7) the proportion 
of market orders. 
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5. Price dynamics 

5.1. Price dynamics in the heavy traffic limit 

Denote by > 0) the (bid) price process corresponding to the Hmit order book process [q^)t>o- 
As explained in Section 2, s" is a piecewise constant stochastic process which 

• increases by one tick at each event (t"'", V^"'") at the ask for which 

• decreases by one tick at each event {t\'^ ^V^'") at the bid for which 

yb.n 

Jn 



Due to the complex dependence structure in the sequence of order durations and sizes, properties 
of the process s" are not easy to study, even in simple models such as those given in Section 2.3. 
The following result shows that the price process converges to a simpler process in the heavy traffic 
limit, which is entirely characterized by hitting times of the two dimensional Markov Q: 

Proposition 1. Under the assumptions of Theorem 2, 

isZt,t>0)"^ S, on (D([0,oo[,R),Mi), where 

•S"* = ^ lQt=o - ^ 1q|=o (25) 

0<s<t 0<s<t 

is a piecewise constant cadlag process which 

• increases by one tick every time the process Q hits the horizontal axis {y = 0} and 

• decreases by one tick every time Q hits the vertical axis {x = 0}. 

Proof. We refer the reader to Whitt (2002) or Whitt (1980) for a description of the Mi topology. 
The price process s" (rescaled in time) can be expressed as 

•5"t = ZJ l*fe(X",Q" ,...,Q" )(r»).(Oa)<0 - l*fc(X",Q" ,...,Q" )(r^).(l,0)<0- 

where r^,Q"n are defined in the proof of Theorem 2. There we showed that 

k 

vfc > 1 , (X" , . . , tj: , , . . . , g';„ ) (x, ^ , . . , r,. , , . . . , q.j . 

As shown in the proof of Theorem 2, {X,Qt^, ...,Qt-^^ ) lies, with probability 1, in the set of 
continuity points of '^k for each A; > 1 so 

(X" , q;;„ ,...,Q^,„j''^^^,iX,Qr„...,Qr,_,). 

Applying Proposition 2 (see Appendix) and the continuous mapping theorem (Billingsley 1968, 
Theorem 5.1) then shows that 

l*fe(X",Q;^„,...,Q;;„ )(t^).(0,1)<0 =^ i*fc(Jf,Qri ,...,Q,^_-^)(rfe).(0,l)<0 

1 fc — 1 

The sequences of processes Z]r"<t l*fc(^".Q"n- -,Q'n )(r").(o,i)<o and 

^1 ■^fe-i 

X]^"<t l*fe(x",Q"„....,Q"„ )(t").(i.o)<o belong to I?^([0, oo[,]R+), the set of increasing cadlag trajec- 

tories. The convergence on the Mi topology for sequences in is equivalent reduces to the 
convergence on a dense subset including zeros. Therefore, 
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and 



On the other hand, since the set of discontinuities of <:t^9k{x,Qri,....Qr )(Tfc).(o,i)<o ^^'^ 



S is thus the difference between the occupation time of the y axis and the occupation time of 
the X axis by the Markov process Q. In particular, this resuh shows that, in a market where order 
arrivals are frequent, distributional properties of the price process s" may be approximated using 
the distributional properties of the limit S. We will now use this result to obtain some analytical 
results on the distribution of durations between price changes and the transition probabilities of 
the price. 

5.2. Duration between price moves 

Starting from an initial order book configuration Qq = {x,y), consider the hitting times 
• the next price increase occurs at the first hitting time of the x-axis by Q: 



which has the same law as the first exit time from the positive orthant of a two-dimensional 
Brownian motion with drift. Using the results of Iyengar (1985), Mctzlcr (2UIU), Zhou (2001) we 
obtain the following result which relates the distribution of this duration to the state of the order 
book and the statistical feature of the order flow process in the case wbalanced order flow where 

W = v^ = o. 

Proposition 2 (Conditional distribution of duration betw^een price changes). In a bal- 
anced order flow where V°- = V^ = ^ the distribution of the duration r until the next price change, 
conditonal on the current state of the bid and ask queues, is given by 




ra = inf{t>0,Q^ = 0} 



• the next price decrease occurs at the flrst hitting time of the x-axis by Q: 



n = mf{t>0, Q^ = 0}. 



The duration r until the next price changes is then given by 



T = TaAn 



nr>t\Q'o = x,Q^, = y] 




(2n + l) 



1 



sm 



(2n-hl)7r6', 



a 



'0 



(/(;.„-i)/2(t:) + -^(^„+i)/2(t:)) 



where ly, 



n 



{2n + l)'K / a. In is the nth Bessel function, 
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IT 

« = < 77 



TT + tan ^ ( 



P 



p>0 

/O = and 

p<0 



vr + tan 



x-py 



vr 
2 



tan ( j 

x-py 



x<py 

x = py (26) 

a; > 



In particular, t is regularly varying with tail index 



vr 

2a' 



Proof. When V°- = V^ = the process Q behaves hke a two-dimensional Brownian motion Z 
with covariance matrix given by (18) up to the first hitting time of the axes, so the distribution of 
the duration r has the same law as the first exit time of Z from the orthant: 

r = inf{t>0,Q;' = or = 0} 

Using the results of Iyengar (1985), corrected by Metzler (2010) for the distribution of the first 
exit time of a two-dimensional Brownian motion from the orthant we obtain the result. 
A result of Spitzer (1958) then shows that 

poo 

¥.[T^\Ql = x,Ql = y]= / t^-'¥[T>t\Ql = x,Ql = y\dt<oo 

TT 

if and only if /3 < 7r/2a, where a is defined in (26). Therefore the tail index of r is — . This result 

does not depend on the initial state {x,y). 

• If p = 0, the two components of the Brownian Motion are independent and r is a regularly- 
varying random variable with tail index 1. This random variable does not have a moment of order 
one. 

TT 

• If /9 < 0, — > 1 and r has a finite moment of order one. In practice, p ~ —0.7; this means that 
if /ia = and /ib = 0, the tail index of r is around 2. 

TT 

• When p > 1, — < 1. The tail of r is very heavy; r does not have a finite moment of order one. 

2a 

For all high frequency data sets examined, the estimates for Pa^P-t are negative (see Section 4.3); 
the durations then have finite moments of all orders. 

Remark 2. Using the results of Zhou (2001) on the first exit time of a two-dimensional Brownian 
motion with drift, one can generalize the above results to the case where {V^,V"-) ^ (0,0): we 
obtain in that case 

2^aixi+a2X2+ait-rQ/2t 



nr>t\Q'o = x,Q^, = y] 
where 9o,a are defined as above, tq 

gn{0) 



at 

U and 



nvrfc 



a 



sm 



rnr 



a 



9n 



(27) 



di = (ai \/%,Va + pa2Vb\/xCj , = (^pa-i \J XaVa + a2Vb\/xCj 

PaVhVb + PbpVaVXa PPaV^Vb + PbVa\/Ki 



ai 



{I - p'^)alXa^/Kvb 



02 



and 



at 



oi + a2 4- 2pai V Xa\bVaa2Vb 



{I - p'^)alXb^/X~aVa 

- aipa- a2Pb- 



(28) 
(29) 

(30) 
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5.3. Probability of a price increase 

A useful quantity for short-term prediction of intraday price moves is the probability y) that 

the price wih increase at the next move given x orders at the bid and y orders at the ask; in our 
setting this is equal to the probability that the ask queue gets depleted before the bid queue. 

In the heavy traffic limit, this quantity may be represented as the probability that the two- 
dimensional process {Qt,t > 0), starting from an initial position {x,y), hits the horizontal axis 
before hitting the vertical axis: 

p-^x, y) = F[Ta < n\{Ql QS) = {x, y)]. 

Since this quantity only involves the process Q up to its first hitting time of the boundary of the 
orthant, it may be equivalently computed by replacing Q by a two-dimensional Brownian motion 
with drift and covariance given by (17)-(18). 

However, when V"- = V'^ = 0, one has a simple analytical solution which only depends on the size 
X of the bid queue, the size y of the ask queue and the correlation p between their increments: 

Theorem 3. Assume V + V'' < 0. Then p"^ : — ^ [0,1] is the unique bounded solution of the 
Dirichlet problem: 

1 ox'^ 1 oy^ oxoy ox ay 

(31) 

with the boundary conditions 

Vj;>0, p''p(x,0) = 1 and Vy > 0, p"P(0,y) = 0. (32) 
When V"- = V^ = 0, Pup{x,y) is given by 



arctan(./f±f ^"" 



1 - - i-P 

2 arctan( 



Pup(x,y) = 2 ' ^^v^,^, ^ (33) 



where Xa,Xb,Va and Vb are defined 1 and 2. 

Proof. Using the results of Yoshida and Miyamoto (1999), the Dirichlet problem (31)-(32) has 
a unique positive bounded solution u G C^(]0,oop,]R+) n C^(]R^,]R+). Application of Ito's formula 
to Mt = u{Q\, Q1) then shows that the process M"^ stopped at r is a martingale, and conditioning 
with respect to {Qq,Qq) = {x,y) gives u{x,y) =p^^{x,y) 

Assume now V"' = V'' = 0. Using a change of variable xt-^ x y/XbVb and y^y V^Va, one only 
needs to consider the case where \/^Vi, = \f\~aVa- 

Up to the first hitting time of the axes, {Q^t > 0) is identical in law to Q = AB where 



A 



cos(/3) sin(/3) 
sin(/3) cos(/3) 



with 13 satisfying p = sin(2/3), /3 < 7r/4 and B a standard planar Brownian Motion with identity 
covariance. Using polar coordinates {x,y) = (rcos0,rsin0) and setting, for r > 0,0 € [0,27r[, 

• ei =^-^1,0) = (-sin(/3),cos(/3)) 

• e2 = ^-i(0,l) = (cos(/3),-sin(/3)) 
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then 

ct>{r,e):=p-^{T ^-^(cos(g),sin(g))=p"^f . (cos(/3 + g), sin(g - 

Vcos''(p) — sm [13) 

is a solution of the Dirichlet problem 

I d d(h 1 d'^d) 

-rd-r^V^7^d$ = ''^ {ir,9),r>0,ee{0,./2)} 

in the cone C = {{r,0), r > 0, g] — /3, | — f3[}, with the boundary conditions 0(]Rei) = 1 and 
(j)(Re2) = 0. A solution, which in this case does not depend on r, is given by 

" ^ (-6' + 7r/2 + arcsin(p)/2), 



7r/2 + arcsin p 



where p is the correlation coefficient between the bid and ask queues. Using the results of Yoshida 
and Miyamoto (1999), this problem has a unique bounded solution so finally 

y) = , ' . f vr/2 + arcsin(p)/2 - arctan(-^^"(^^^*""(^/^) " 



7r/2 + arcsin /3 \ cos(/3 + arctan(?//x)) 

Remark 3. When \/^Va = \/^Vb, the probability p^P{x,y) only depends on the ration y/x and 
on the correlation p 

^ arctan(yi±^^) 

Pup {x,y) = - , (34) 

^ 2arctan(Jg) 

and when /) = (which is the case for some empirical examples, see Section 4.3), 

2 y 
p^^(x,y) = — arctan( — ). 

TT X 

Figure 15 displays the dependence of the uptick probability on the bid-ask imbalance variable 
9 = arctan(y/x) for different values of p. 
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< arctan(y/x) < 7i/2 



Figure 15 p^^ as a function of tlie bid-ask imbalance variable 6 — arctan(j//2:) for p 
(green line) and p = —0.9 (red line). 



= (blue line), p= -0.7 
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6. Appendix: Technical Proofs 

6.1. A Ji-continuity property 

Lemma 2. Let t : Z)([0,oo),]R^) i— )■ [0, cxd[ be the first exit time from the positive orthant. The map 

G: (L'([0,oo),M'), Ji) ^ M (35) 

W — )■ 1w(t(w)).(0,1)<0- 

is continuous on the set {uj G C([0, oo),M^\{(0, 0)}), t{uj) < oo}. 

When t{uj) < oo, G{u}) = 1 indicates that w first exits the orthant by crossing the x-axis. To prove 
this property, first note that 

C([0,oo),M^{(0,0)})= Uc^o([0,oo),M^i?(0,l/n)). 

n>l 

Let Wo G C([0,oo),M2\{(0,0)}). There exists n e N such that ojq i B{0,l/n). Let e > such that 
e + ^wo(^) +^woo>'(^) ^ ^/^^ where tj^^ is the modulus of continuity of uj. Let uj' £ -D([0,oo),M^) with 
dj-^{uJo,uj') < e. There exists A : [0,T] — )• [0,T] increasing such that: 

1 o A — a;| loo < e and ||A — e||oo- 

Without loss of generality, one can also assume, by continuity of r on {D,Ji), that 

\t{ujo) -t{uj)\ < e. 

Now, we will show that |a;o(T(c<^o)) - w'(r(a;'))| < e + r/^o(e) + ?7„(,oA(e): 

\ujo{t{ujq))-u}'{t{uj'))\ = \uo{t{uo))-ujooX{t{u}'))+ujooX{t{uj'))-ujooX{t{ujo)) + ujooX{t{u}o))-uj'{t{uj'))\, 
therefore 

\ujo{t{ujo)) - uj' {t{uj'))\ < \ \uJao X-uj' Woo + \ujoo X{t{uj')) - ujq o A(r(wo))| + \^o° Ht{ujo)) - ujo{t{ujo))\ 

Since e + r?„o(e) + ??^ooA(e) < 1/"- and ujq ^ B{0, l,n), lr(c^o) (o,i)<o = lr(a;').(o,i)<o, which completes 
the proof of the continuity of the map G on the space C([0,oo),]R^\{(0,0)}). 

6.2. Continuity of ^: proof of Theorem 1 

To study the continuity of the map we endow Z?([0,oo),M^) with Skorokhod's Ji topology (see 
Lindvall (1973), Whitt (1980)). Let Ar the set of continuous, increasing functions A : [0,r] [0,T] 
and e the identical function on [0,T]. Recal that the following metric 

{uji,uj2) = inf (||a;2 A-a;i||oo + ||A-e||oo) • 

defined for uji,uj2 G Z)([0,r],M^), induces the Ji topology on D{[0,T],R'^), and uj,, uj in 
(L'([0,oo),M2), Ji) if for every continuity point T ofuj, uj^^uj in {D{[0,T],R'^), J^). 

The set (M^) is endowed with the topology induced by 'cylindrical' semi- norms, defined as 
follows: for a sequence (i?")„>i in (M^) 
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D([0, cxd),]R^) X (M^) X (M^) is then endowed with the corresponding product topology. The goal 
of this section is to characterize the continuity set of the map 

introduced in Definition 1. Let us introduce C([0, oo),]R^\{(0, 0)}) be the space of continuous planar 
paths avoiding the origin: 

C([0,oo),r\{(0,0)})=UCo([0,r],r\S(0,l/n)). 

ri>l 

Lemma 3. Let cu e C{[0,oo),R'^\{{0,0)}). Then the map 
*i :L'([0,oo),]R^) XM+ xM+ ^ i:»([0, oo), M^) (36) 

{uj,Ri,Ri) l-> CJ + l[r(w),oo) (^1ct6(w)=t(w)(-Ri —^t{uj}) + IfTa (w)=t(w) (-Rl — ^t{uj})^ , 

where 

(7fc(a;)=inf{t>0, wt.(0,l)<0}, o-a(w) = inf{t > 0, a;t.(l, 0) < 0} and t{uj) = at{uj) Aaa{uj). 

is continuous at oo with respect to the following distance on (£)([0, cxo],]R^) x x M+)." 

d{{uj,Ri,Ri),{uj\R[,R!^)) = dj^{Lo,J) + \Ri- R:^\ + \Ri- Ri\ and: 

Proof. Let {u}o,Ri,Ri) G C([0,oo),M2\{(o,0)}) x M^, {u}',R[,R[) E D{0,m'^) x M^. Since u>o € 
C([0,oo),M^\{(0,0)}), there exists n > such that ojq ^ B{0, 1/n). Let < e < 1/n such that 

d{iu;o,Ri,Ri),iuj',R[,R[))<e. 

Since dJ-^^{uJo,uJ') < e, there exists A : [0,T] — [0,T], non-decreasing such that: 

— e||oo<e, and | Icjq o A — a;||oo < e- 

By continuity of r for the Ji topology Whitt (2002) [Theorem 13.6.4], one can also assume, 
without loss of generality, that 

\t{ujo o A) — t{uj')\ < e. 



Moreover, since the graph of Uq does not intersect with i?(0,l/n) and e < 1/n, lT{uio)=aa{t^o) ~ 
'i-T{u>')=c7aiLj')- Now define A' by 

A^[0,T] ^ [0,r] (37) 
t(wooA) 



Then 



t(uj) e 



r(wooA) t(wooA) 
On the other hand 

IIwqo A' - w||oo = IIwqo A' - A + a;oo A - wiloo < 11^0° A' - A||oo + e < r]^o°xie) + e, 
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where ij^^ox is the modulus of continuity modulus of ojq o A. Therefore, sincelT-((^uoAE) = It(w') by 
definition of A' and 

* 1 ( Wo , fli , ^1 ) o A' - ^' 1 (a;' , , ^'i ) ) = Wo o A' - 

Thus A*^ satisfies 1 1 A' - ell <e( ^'"''+^ ) and 

1 1 *i (wo , «i , i?i ) o A^ - ^-1 (w' , , ) ) 1 1 oo < r?^oo A (e) + e + 2e 

which proves that {uq, Ri, Ri) is a continuity point for 'I'l. 
For k>2, define recursively the maps 

'ffc:L>([0,oo),M') xM^ xM^ ^ £'([0, oo), M^) (38) 

{UJ,{R„R^),>1) ^ ^l{^k-li^^,{R^,R^)^=l..k-l),Rk,Rk). 

To simplify notation we will denote the argument of ^'^ as {uj, R, R){= (to, {Ri, Ri)i>i) although it 
is easily observed from (38) that ^'fe only depends on the first k elements {Ri,Ri)i=i,,k) of R,R. 

Lemma 4. If {uj,R,R) e C{[0,oo),M.'^\{{0,0)}) xWl xR^ such that 

{0,0)^^,{oj,R,R){[0,oo) ) (39) 

then is continuous at (uj,R,R). 

Proof. Let {Ri, Ri)i>i, (i?-, i?-)i>i, two sequences of random variables on and define 

n^iR, R) = n^-^o^. (<^([0, oo),M^{(0, 0)}), i?, R) 

where we have set = Id. Consider ooq E Qk{R,R), and lu S Z?([0,T],M^), such that: 

dj^{uo,u;)+ sup \Ri-Rl\+ sup \Ri-R[\<e. 

i—l..k i—l..k 

An application of the triangle inequality yields 

dJ^{^k{0J0,{R^,R^)),^k{0J',{R',,Rfi})) 
< dj^ i^kiooo, {R^,R^)),^k{uJ', (i?„ R,))) + dj, k{oj' , {R, , R,)) , ^ k{uj' , (iJ^, R!^)) 

where the last term converges to zero when e goes to zero by continuity of ^'i. 
We can now prove Theorem 1. 

Proof of Theorem 1. Since uj is continuous, the jumps of ^'(w, R, R) correspond to the first exit 
times from the orthant of the paths ^'^(w, i?, i?). Therefore, if (-R„)„>i, (-Rn)n>i have no accumula- 
tion points on the axes, the paths ^{uj, R, R) only has a finite number of discontinuities on [0, T] for 
any T > 0. So, for any T > 0, there exists k{T) such that ^ = ^k{T) - Then thanks to Lemma 4, ^ is 
continuous on the set of continuous trajectories whose image has a finite number of discontinuities 
and does not contain the origin. 
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6.3. Functional central limit theorem for the net order flow 

Proposition 3. Let {T^'" ,T^'")iyi and {V°"" ,V^'''")i>i be stationary arrays of random variables 
which satisfy Assumptions 1 and 2. Let (N^'",t > 0) and {Nf'^\t > 0) be the counting processes 
defined in (15). Then 

nt -rra,n nt -rrb^n \ 

E^'E^ ^ (^B,+t{X'^W,\''V^)) (40) 
^-^ Jn ^-^ Jn \ n^cx> \ / f>o 

i—1 ^ i—1 ^ I 

/ t>0 

where B is a standard planar Brownian motion and 

^^-[pVX^^v^v, \%l J' ^^^^ 
Proof: First we will prove that the sequence of processes 

4 (sBt + t(A'^y-,AV^)) 

\/n ^-^ \/n I n-i-oo V /t 



Jn — ' wn / n-i-oo \ / t>o 

1=1 * i=i ^ I - 

/ i>0 

weakly converges in the Ji topology. Using the Cramer- Wold device, it is sufficient to prove that 
for (a,/3) gM2^ 



a 



I- "J T/;o,n I- -J yo.n \ , 

V^ + /3V^ {a\^V- + p\''V''))t + J{a^X-vl + P^\Hl + 2paPvaVbV¥\^)Bt 

— ^ \/TI — ^ \/TI ' n.^on 
1 ^ 2=1 ^ 



n—>oo 

t>0 



If A'^ G Q and E Q, it is possible to find A such that A"/A E N and AVA E N. Let for all (z, n) E W 



then for all i > 0, 



1^°''^] ya.n [^-^A yb.n [^*] 



i—1 i—1 i—1 



For all n > 0, {W",i > 1) is a sequence of stationary random variables. Therefore, thanks to theorem 
(Jacod and Shiryaev 2003, Chap.VIII, Thm 2.29, p.426), and the fact that 

OO 

Y&riW^) + 2 ^ cov{Wl\ WJ') ™ cj2, (42) 

1=2 

S!=T' t^O] converges weakly to a Brownian motion with 

V ^ / n>l 

volatility \/Acr. If (A^ A**) ^ there exists (A^, A^J„>i such that 

X^X'^eQ and |A::-A1<-, |A:;-A''|<i 

n n 



A A" 

As above, one can define an integer A„ such that € Q and e Q. Let for all ii.n) G N' 
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One has for all t>0, 

ttVn {rtVn tt U 

Moreover 

therefore the convergence above holds even if or A;, are not rationals. On one hand, 

var(W,") = var + ... + ,J + /3(F^^5/,„,,_,hi + + - + ^S/aJ) 

+ 2a/3cov + ^aIva^- ^(A^/A„)(i-i)+i- + K^/xS) ■ 



On the other hand, for alH > 2, 



co^m, wn = a^cov (vr + ... + v,^:,,^, , v- + - + 

+ /3C0V Vi' +--- + V,\,. ^,V^,6/^ w n^i + ••• + -/i 

+ a/5cov (f-" + ... + ^rA^/A„), + - + ^a^/aJ 

+ a/3cov (f,^- + ... + V^(^4"/,„,, V(a':/a„),-i).i + - + n",r/A„) ■ 

Therefore 

var(Wi") + 2 V cov(Wi", W,") = var(yi"'") ^ + 2 V cov(yi"'", F,"'") ^ 

— ^ — ^ Ari 

+ var(yi'''") ^ + 2 V cov(Fi'''", F/'") ^ 

A„ ^" 



+ 



2a/3cov (f-... + F-, y^... + J 

i/A„)(«-l) + l ^ ••• ^ "^Afcz/Any 

+ 2a/35:cov (y^"... + ^4"/,^, ^(a':/a„),-i)+i + - + V^^,,^) 

A simple calculation shows that 
2a^C0Y[y,^ ••• + ^AS/A„'^i' ••• + ^Ai/Aj 

+ y: -V {vr... + ^a1m„ , ^ir/A„),-i).. + - + nS/A J 



+ 2a/35:cov(K^-... + y^;;,„,i^^r/. .,,+... + v^'^ 

i=2 



i=2 
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1—2 



Therefore 

lim var(t^/^) + 2ycov(T^i",Wr) = «^^a + /3^^' + 2^ 

1.1— ^no ^ — ^ A A \ 



1=2 

where p is given in (11) and 



lim E[Wr] = a^y^ + /3^y^ 



which completes the proof of the convergence in (42). The law of large numbers for renewal processes 
implies that the following sequence of processes converges to zero in the Ji topology Iglchart and 
Whitt (1971): 

(iV„T)*>o™([A'^i])*>o, and (iV^i"),>o "^°°([A4]),>o, 

nt -^a.n nt ^o.n \ 

^ ^0 in the Ji topology. 



t>0 



6.4. Identification of the heavy traffic limit Q 

Lemma 5. The process Q is a Markov process with values in — {(0,0)} and infinitesimal gen- 
erator given by (13)- (14) and domain 

dom(g) = {/iGC2(]0,oo[x]0,oo[,M)nC"(M^,M), Vx > 0, Vy > 0, 

hix,0)= I h{g{ix,0),{u,v)))F{du,dv)=0, h{0,y) = J^, h{g{iO,y), {u,v)))Fidu,dv) = 0} 

Proof. We use the explicit construction of Q from the planar Brownian motion X using the 
maps defined in the proof of Theorem 2. First, let us show that the process Q is a Markov 
process. Let < ti < ... < tn < t, and k such that Tk <tn< Tk+i. Let A a measurable set. The 
following equations: 

nQt G ^IQt,, QtJ = P[Qt G ^IQt,, gt„, K < t < rfc+i}]P[t„ < t < Tk+i\Qt,....MtJ\ 

+ ¥[Q, e A|Q,,, Q,„,{t > n+i}nt > n+i\Qt, , QtJ 
= P[Qt £ A\Qt„,{t„ < t < n+i}]F[t„ < t < Tk+i\QtJ 
+ P[Qt G A\Qt„,{t > n+,}nt > n+i\QtJ 
= F[Q,gA\Q,J 

prove that the process Q is a Markov process. Let us now compute its infinitesimal generator Q. 
The domain dom(^) of G consists of all functions h S C^(M^) verifying for all {x,y) G 

y nHQt)-h{Q,)\Q, = {x,y)] ^ 

lim < oo. 

t->o t 

For X > 0, and y > 0, a classical computation shows that if /i € C^(]R^), 

nh{Q,)\Q, = {x,y)] = h{x,y) + t (^A.^^- + X.V- + ^q^2 + ^Qy2 +^P^'>^^^>'q^) 
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which leads to equation (13). On the other hand, for all (x,y) € M^, 

IE[/i(Qt)|Qo = (a;,0)]= / ¥.[h{Qt)\Qo+ = g{{x,^),{u,v))]F{du,dv) 

(E[/i(Q,) IQ0+ = g{{x: 0), (n, v))] - h{g{{x, 0), (n, v))) + h{g{{x, 0), {u, v)))) F{du, dv) 
{tgh{g{{x, 0), {u, v))) + h{g{{x, 0), (n, «)))) F(d^x, dv) + o{t). 



nHQt)\Q, = (0,y)] = E[/i(Q,)|Qo+ = ff((0,y), 

(E[/i(Q,)IQo+ = 5((0,y), Kt-))] - h{g{{0,y), {u,v))) + /i(ff((0,y), (n, t;)))) dt;) 
{tghigiiO, y), {u, v))) + h{g{{0, y), (n, v)))) F{du, dv) + o(i). 



As t — )• 0, these equations lead to 

nh{Qt)\Qo = {x,0)]-h{x,0) 



gh{g{{x,0),{u,v)))F{du,dv) 



+ - [ h{g{{x,0),{u,v)))-h{x,0)F{du,dv) + o{l). 



E[MQt)IQo = (o,i/)]-Mo,y) 
t 



gh{gi{0,y),{u,v)))F{du,dv) 



+ 7 / {h{g{{0,y),{u,v)))-g{0,y))F{du,dv) + o{l). 

Therefore, the domain dom(^) of g consists of all function h G C^(]R^ x M^) such that 
/i(x,0)= / hig{{x,0),iu,v)))Fidu,dv) = 0, /i(0,y)= / h{g{{0,y),iu,v)))F{du,dv) = 0, 

When these 'boundary conditions' are verified, the above limits exist and we have 

gh{x,0)= / gh{g{{x,0),{u,v)))F{du,dv), 

gh{0,y)= I ghig{iO,y),{u,v)))Fidu,dv), 

which concludes the proof. 
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